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Abstract 

We consider M = 1 supersymmetric QCD with the gauge group U(iV) and 
Nf = N quark flavors. To get rid of flat directions we add a meson superfield. 
The theory has no adjoint fields and, therefore, no 't Hooft-Polyakov monopoles 
in the quasiclassical limit. We observe a non-Abelian Meissner effect: conden- 
sation of color charges (squarks) gives rise to confined monopoles. The very fact 
of their existence in J\f = 1 supersymmetric QCD without adjoint scalars was 
not known previously. Our analysis is analytic and is based on the fact that the 
]\[ = 1 theory under consideration can be obtained starting from N = 2 SQCD 
in which the 't Hooft-Polyakov monopoles do exist, through a certain limiting 
procedure allowing us to track the status of these monopoles at various stages. 
Monopoles are confined by BPS non-Abelian strings (flux tubes). Dynamics of 
string orientational zero modes are described by supersymmetric CP{N — 1) 
sigma model on the string world sheet. If a dual of J\f = 1 SQCD with the 
gauge group U(iV) and Nf = N quark flavors could be identified, in this dual 
theory our demonstration would be equivalent to the proof of the non-Abelian 
dual Meissner effect. 



1 Introduction 



Seiberg and Witten demonstrated [1] that the dual Meissner effect takes place in Af = 
2 Yang-Mills theories. Shortly after [1], Hanany, Strassler and Zaffaroni discussed 
[2] formation and structure of the chromoelectric flux tubes in the Seiberg-Witten 
solution. Their analysis showed that details of the Seiberg-Witten confinement are 
quite different from those we expect in QCD-like theories. The confining strings in the 
Seiberg-Witten solution are, in fact, Abelian strings of the Abrikosov-Nielsen-Olesen 
type [3]. The "hadronic" spectrum in the Seiberg-Witten model is much richer than 
that in QCD (for a review see e.g. [4].) The discovery of non- Abelian strings [5, 6] and 
non- Abelian confined monopoles [7, 8] was a significant step towards QCD. They were 
originally found in Af = 2 models which are quite distant relatives of QCD. To get 
closer to QCD one needs to have less supersymmetry. Another conspicuous feature 
of Af = 2 Yang-Mills theories which drastically distinguishes them from QCD-like 
theories is the presence of scalar and spinor fields in the adjoint representation. 

To advance along these lines it is highly desirable to break Af = 2 down to Af = 
1 and get rid of the adjoint superfield by making it very heavy, without destroying 
non- Abelian strings and confined monopoles. A partial success in this direction was 
reported in Ref. [9]. Adding a mass term to the adjoint superfield of the type 5W = 
\xA 2 breaks Af = 2 . As long as the mass parameter \i is kept finite, the non- 
Abelian string in this Af = 1 model is well-defined and supports confined monopoles. 
However, at /i — > oo, as the adjoint superfield becomes heavy and we approach the 
limit of Af = 1 SQCD, an infrared problem develops. This is due to the fact that 
in Af = 1 SQCD defined in a standard way the vacuum manifold is not an isolated 
point; rather, there exists a flat direction (a Higgs branch). On the Higgs branch 
there are no finite-size BPS strings [10]. Thus one arrives at a dilemma: either one 
has to abandon the attempt to decouple 1 the adjoint superfield, or, if this decoupling 
is performed, confining non- Abelian strings cease to exist [9]. 

In this paper we report that a relatively insignificant modification of the bench- 
mark Af = 2 model solves the problem. All we have to do is to add a neutral meson 
superfield M coupled to the quark superfields through a superpotential term. Act- 
ing together with the mass term of the adjoint superfield, it breaks Af = 2 down 
to Af = 1 . The limit \i — > oo in which the adjoint superfield completely decouples, 

1 Below we use the word decouple in two opposite meanings: first, if a field becomes very heavy 
and disappears from the physical spectrum, so that it can be integrated out; second, if all coupling 
constants of a certain field vanish so that it becomes sterile. With regards to the adjoint fields 
decoupling means making them very heavy. With regards to the meson superfield M decoupling 
means sterility. Each time it is perfectly clear from the context what is meant. We hope this will 
cause no confusion. 
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becomes well-defined. No flat directions emerge. The limiting theory is Af = 1 SQCD 
supplemented by the meson superfield. We show that it supports non-Abelian strings. 
The junctions of these strings present confined monopoles, or, better to say, what be- 
comes of the monopoles in the theory where there are no adjoint scalar fields. There 
is a continuous path following which one can trace the evolution in its entirety: from 
the 't Hooft-Polyakov monopoles which do not exist without the adjoint scalars to 
the confined monopoles in the adjoint-free environment. As far as we know, this is 
the first demonstration (in fully controllable weak coupling regime) of the Meissner 
effect in Af = 1 theories without adjoint superfields. If a dual of Af = 1 SQCD with 
the additional meson superfield could be found, in this dual theory our demonstration 
would be equivalent to the proof of the non-Abelian dual Meissner effect. 

The organization of the paper is as follows. In Sect. 3 we review the benchmarks 
A" = 2 super- Yang-Mills theory with the gauge group U(iV) and Nf = N quark 
flavors. We introduce the Fayet-Iliopoulos (FI) term [11] in the U(l) subgroup, 
crucial for the string construction, and a meson superfield M, coupled to the quark 
superfields through a cubic superpotential. We add the mass terms to the adjoint 
superfields. The latter two terms in the superpotential break Af = 2 . In Sect. 3 we 
discuss the spectrum of elementary excitations, in particular, in the limit // — > oo. 
We show that the limiting theory is essentially Af = 1 SQCD. The only distinction is 
the meson superfield which survives in the limit /J, — > oo. The vacuum of this theory, 
which will be referred to as M model, is isolated (i.e. there are no flat directions). As 
usual, we construct non-Abelian strings and determine the world-sheet theory. This 
is the contents of Sect. 4. One of the crucial points of our analysis is determination 
of the fermion zero modes. To count these modes we engineer an appropriate index 
theorem (Sect. 5). This theorem applies to the two-dimensional Dirac operator which 
we encounter in the string analysis. In Sect. 5.1 we derive the index v = 4N. We 
observe four supertranslational zero modes and 4(JV — 1) superorientational modes. In 
Sect. 6 we discuss how the monopoles evolve when we vary the adjustable parameters 
of the M model: from the 't Hooft-Polyakov limit to the limit of confined monopoles 
in highly quantum regime in Af = 1 SQCD. In Sect. 7 the same issue is discussed 
from the brane perspective. Section 8 summarizes our findings. Finally, in Appendix 
we present explicit expressions for the fermion zero modes in the case of two flavors. 

2 From AT = 2 SQCD to AT = 1 

To begin with, let us briefly review Af = 2 supersymmetric QCD. The gauge symme- 
try of our benchmark model is SU(iV) xU(l). It has Nf = N matter hypermultiplets. 
The field content of this model is as follows. The Af = 2 vector multiplet consists of 
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the U(l) gauge fields and SU(iV) gauge field A a , (here a — 1, A^ 2 — 1), their 
Weyl fermion superpartners (A*, A 2 ) and (A* a , A 2 "), and complex scalar fields a and 
a a , the latter in the adjoint of SU(iV). The spinorial index of A's runs over a — 1, 2. In 
this sector the global SU(2)r symmetry inherent to M = 2 model at hand manifests 
itself through rotations A 1 <-> A 2 . 

The quark multiplets of the SU(iV)xU(l) theory consist of the complex scalar 
fields q kA and qAk (squarks) and the Weyl fermions ip kA and ipAk, all in the fundamen- 
tal representation of the SU(iV) gauge group. Here k = 1,...,N is the color index 
while A is the flavor index, A — 1, ...,N. Note that the scalars q kA and q kA form a 
doublet under the action of the global SU(2) R group. 

In addition, we introduce the Fayet-Iliopoulos D-term for the U(l) gauge field 
which triggers the squark condensation. 

The undeformed Af = 2 theory we start from has a superpotential, 

nV =2 = v^Tr ^-QAQ + QA a T a Q^ + TrmQQ (2.1) 

where A a and A are chiral superfields, the M = 2 superpartners of the gauge bosons 
of SU(iV) and U(l), respectively, while T a are generators of SU(iV) normalized by 
the condition 

Tr r j ]ar j ] b — — 
~ 2 

Moreover, m is the quark mass matrix, a numerical N x N matrix m\ (to be elevated 
to superfield matrix later on). We write the quark superfields Q kA as N x N matrices 
in color and flavor indices. The trace in (2.1) runs over the appropriate indices. 

Now we deform this theory in two ways each of which breaks M = 2 supersym- 
metry down to N = 1 . First, we add superpotential mass terms for the adjoint chiral 
superfields from the U(l) and SU(iV) sectors, respectively, 

^/fy^ + f^ 2 ' (2 - 2) 

where H\ and H2 are mass parameters. Clearly, the mass term (2.2) splits the gauge 
M = 2 supermultiplets, breaking M = 2 supersymmetry down to M = 1 . As will 
be discussed later in detail, in the large-// limit the adjoint multiplets decouple and 
then we recover M = 1 SQCD with Nf = N flavors. This theory has a Higgs branch 
(see, for example, [12]). The presence of quark massless states in the bulk associated 
with this Higgs branch obscure physics of non-Abelian strings in this theory [9]. In 
particular, the strings become infinitely thick. 

Can one avoid this shortcoming? The answer is yes. To this end we introduce 
another J\f = 2 breaking deformation. Namely, we uplift the quark mass matrix 
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to the superfield status, 



where M represents N 2 chiral superfields of the mesonic type (they are color- singlets). 
With this uplifting we have to add a kinetic term for Aff , 

^Mkin = Jd 4 xd 2 9d 2 9 — Tr MM , (2.3) 

where h is a new coupling constant (it supplements the set of the gauge couplings). 
At h = the matrix field M becomes sterile, it is frozen and in essence returns to 
the status of a constant numerical matrix as in Ref. [9]. The theory acquires flat 
directions (a moduli space). With nonvanishing h these flat directions are lifted and 
M is determined by the minimum of the scalar potential, see below. 

The elevation of the quark mass matrix to superfield is a crucial step which 
allows us to lift the Higgs branch which would develop in this theory in the large \x 
limit if M were a constant matrix. 

The bosonic part of our SU(iV)xU(l) theory has the form 



+ Tr \V„q\' + Tr |V„?|" + ^ 



d„M 



+ -\d fl M a \ z + V(q,q,a a ,a,M iJ ,M a ) 



(2.4) 

Here is the covariant derivative in the adjoint representation of SU(2), while 

V» = d»- l -A,-iAlT a . (2.5) 



Moreover, the matrix M^ can be always decomposed as 

M^ = -5^ M° + (T a )i M a . (2.6) 
2 

We use this decomposition in Eq. (2.4). The coupling constants gi and g 2 correspond 
to the U(l) and SU(A^) sectors, respectively. With our conventions the U(l) charges 
of the fundamental matter fields are ±1/2. 

The potential V(q A , q^, a a , a, M°, M a ) in the Lagrangian (2.4) is a sum of various 
D and F terms, 



V(q A ,q A ,a a ,a,M°,M a ) = % ( \ f abc a b a c + Tr qT a q - Tr qT a ~q 



2 U 



4 



+ IcTrfg-Tr^-A^ + f 



+ i 



Tr qq + v iV/xia 



+ 2 Tr 



2TrgT a g + ^2/z 2 a a 



> + 2 T a a a )g + -^q(M° + 2T a M a ) 
v2 



+ 



(a + 2 T a a a )q + -^q(M° + 2T a M a ) 
v2 



- |Tr ggT + h\TiqT a q 



~|2 



(2.7) 



where / abc stand for the SU(iV) structure constants. The first and second terms 
here represent D terms, the next two terms are F4 terms, while the term in the curly 
brackets represents the squark F terms. Two last terms are F terms of the M field. In 
Eq. (2.7) we also introduced the FI D-term for the U(l) field, with the FI parameter 
£. Note that the FI term does not break M = 2 supersymmetry [2, 13]. The three 
parameters which do break Af = 2 down to M = 1 are fj,i, fi 2 and h. 

The FI term triggers the spontaneous breaking of the gauge symmetry. The 
vacuum expectation values (VEV's) of the squark fields can be chosen as 



(q kA ) = ^ 



( 1 



V 








(t A ) = 0, 



k = 1,...,N, A=1,...,N, 
up to gauge rotations. The VEV's of the adjoint fields vanish, 

(a a >=0, (a) = 0, 

and so do those of the M fields, 

(M a )=0, (M°)=0. 



(2.8) 



(2.9) 



(2.10) 



The color-flavor locked form of the quark VEV's in Eq. (2.8) and the absence of 
VEV's of the adjoint scalar a a and the meson scalar M a in Eqs. (2.9), (2.10) results 
in the fact that, while the theory is fully Higgsed, a diagonal SU(iV)c + .F symmetry 
survives as a global symmetry. Namely, the global rotation 

q^UqU' 1 , a a T a ^Ua a T a U-\ M^U~ X MU, (2.11) 

where U is a matrix from SU(iV) is not broken by the VEV's (2.8), (2.9) and (2.10). 
This is a particular case of the Bardakci-Halpern mechanism [14]. The presence of 
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this symmetry leads to the emergence of orientational zero modes [6] of the Z N strings 
in the model (2.4). 

Note that the vacuum expectation values (2.8), (2.9) and (2.10) do not depend 
on the supersymmetry breaking parameters /ii and H2- This is because our choice of 
parameters in (2.4) ensures vanishing of the adjoint VEV's, see (2.9). In particular, 
we have the same pattern of symmetry breaking all the way up to very large values 
Hi and /i 2 , where the adjoint fields decouple. 

With N matter hypermultiplets, the SU(iV) part of the gauge group is asymp- 
totically free, implying generation of a dynamical scale A. In the ultraviolet (UV) 
we start with a small g%, and let the theory evolve in the infrared. If the descent to 
A were uninterrupted, the gauge coupling g\ would explode at this scale. Moreover, 
strong coupling effects in the SU(A^) subsector at the scale A would break the SU(Af) 
subgroup through the Seiberg-Witten mechanism [1]. Since we want to stay at weak 
coupling, we assume that a/£ ^> A, so that the SU(AQ coupling running is frozen by 
the squark condensation at a small value, namely, 

^ = ln^ + ...»l. (2.12) 
Ngl A 

Now let us discuss the elementary excitation spectrum in the theory (2.4). Since 
both U(l) and SU(AT) gauge groups are broken by the squark condensation, all gauge 
bosons become massive. From (2.4) we get for the U(l) gauge boson mass (we refer 
to it as photon) 




™ P h = #iyy£, ( 2 - 13 ) 
while (A^ 2 — 1) gauge bosons of the SU(AQ group acquire a common mass 

m w = g 2 Ji. (2.14) 

This is typical of the Bardakci-Halpern mechanism. To get the masses of the scalar 
bosons we expand the potential (2.7) near the vacuum (2.8), (2.9), (2.10) and diago- 
nalize the corresponding mass matrix. The A^ 2 components of the 2 A^ 2 -component 2 
scalar q k are eaten by the Higgs mechanism for U(l) and SU(AQ gauge groups. An- 
other A^ 2 components are split as follows: one component acquires the mass (2.13). 
It becomes a scalar component of a massive M = 1 vector U(l) gauge multiplet. The 
remaining A^ 2 — 1 components acquire masses (2.14) and become scalar superpartners 
of the SU(Af) gauge boson in the Af = 1 massive gauge supermultiplet. 



2 We mean here real components. 



6 



Moreover, 6 N 2 real scalar components of the fields q^k, a a , a, M a and M° 
produce the following states: six states have masses determined by the roots of the 
cubic equation 

\l - A 2 (2 + u 2 + 2 7l ) + A,(l + 2 7l + 7- + 2 W ) - 7i ^ = , (2.15) 
for i — 1. Namely, these states form degenerate pairs with the masses 



iV 

"iu(i) =9iU-2&i- (2-16) 

Each root of Eq. (2.15) for % — 1 determines masses of two degenerate states. 

Above we introduced M = 2 supersymmetry breaking parameters cu and 7 asso- 
ciated with the U(l) and SU(iV) gauge groups, respectively, 

7T 2 ~7T (2 ' 17) 

and 

h h , 

i^W (2 ' 18) 

Now we are left with 6 (iV 2 — 1) states. They acquire masses 

rasu(iv) = 92\f£\i , (2.19) 

where each root of Eq. (2.15) for % = 2 determines masses of 2 (A^ 2 — 1) degenerate 
states. 

When the supersymmetry breaking parameters uj{ and 7^ vanish, two mass eigen- 
values (2.16) coincide with the U(l) gauge boson mass (2.13). The corresponding 
states form the bosonic part of the M = 2 long massive U(l) vector supermultiplet 
[13]. The one remaining eigenvalue in (2.16) becomes zero. It corresponds to the 
massless field M° which decouples (becomes sterile) in this limit. With nonvanishing 
values of lo\ and 71 this supermultiplet splits into the massive Af = 1 vector multiplet, 
with mass (2.13), plus three chiral multiplets with masses given by Eq. (2.16). 

The same happens with the states with masses (2.19). If cu's and 7's vanish they 
combine into the bosonic parts of (iV 2 — 1) M = 2 massive vector supermultiplets, 
with mass (2.14), plus the massless field M a . If u;'s and 7's do not vanish, these 
multiplets split into (iV 2 — 1) M = 1 vector multiplets (for the SU(iV) group), with 
mass (2.14), and 3 (iV 2 — 1) chiral multiplets, with masses (2.19). 
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3 Af = 1 SQCD with the mesonic M field 



Now let us take a closer look at the spectrum obtained above, assuming the limit of 
very large Af = 2 supersymmetry breaking parameters 

OJi > 1. 

In this limit the largest masses tou(i) and msu(jv) become 

(largest) 

m U(l) = m U(l) = 

(largest) 2 /o 1 \ 

m SU(W) = m SU(2) ^2 = # 2 ^2 • (3-1) 

Clearly, in the limit /ij — > oo these are the masses of the heavy adjoint scalars a and 
a a . At cjj ^> 1 these fields leave the physical spectrum; they can be integrated out. 

The low-energy bulk theory in this limit contains massive gauge Af = 1 multiplets 
and chiral multiplets with two lower masses mu(i) and two lower masses msu(TV)- 
Equation (2.15) gives for these masses 




(1 ) _ JWVJ 



'U(l) 

for the U(l) sector and 



(2) hNi 



[ l -2^V^Tl) + •••}, (3.2) 



4V) = /f { 1 + ^V / 7 2 (72 + !) + •••}, 

l-^-V^ + T) + •••}, (3.3) 



,(2) 




m SU(V) 

for the SU(iV) sector. 

It is worth emphasizing again that there are no massless states in the bulk theory. 
As we have already mentioned, at h — the theory (2.4) develops a Higgs branch in 
the large-// limit (see, for example, [9]). If h ^ 0, M becomes a fully dynamical field, 
and the Higgs branch is lifted, as follows from Eqs. (3.2) and (3.3). 

At large \i one can readily integrate out the adjoint fields A a and A. Instead of 
the superpotential terms (2.1) and (2.2) we get 

(TtQQ) 2 (TrQT a Q) 2 „ 

W = -^- J — + TtMQQ. 3.4 

4/xi fi 2 
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At /j,i t 2 — * oo the first two terms disappear, we are left with W = Tr M Q Q, and our 
model (2.4) reduces to Af = 1 SQCD with the extra mesonic M field. The bosonic 
part of the action takes the form 

s = I dtx ^|( f ;-) 2+ 4^( f, '"') 2+Tr i v "''i 2+Tr i v ^i 2 



+ h 



^M 



+ - |9 M M a r + ^ (TrgT a g - Tr qT a q)' 

ft Zi 



a 2 

+ ^ (Trgg - Tr gg - JVf) 2 + Tr|gM| 2 + Tr|gM| 2 
8 

+ ||Trgg| 2 + /i|TrgT a g| 2 | . (3.5) 

The vacuum of this theory is given by Eqs. (2.8) and (2.10). The mass spectrum 
of elementary excitations over this vacuum consists of the M = 1 gauge multiplets 
for the U(l) and SU(iV) sectors with masses given by Eqs. (2.13) and (2.14), and 
the chiral multiplets of the U(l) and SU(iV) sectors with masses given by the leading 
terms in Eqs. (3.2) and (3.3). The scale of the theory (3.5) is determined by the scale 
of the theory (2.4) in the M = 2 limit by the relation 

A^ =1 = tfA» . (3.6) 

In order to keep the theory (3.5) at weak coupling we assume that 

( ?2V / e>A A r=i. (3.7) 

Our Af = 1 SQCD with the M field, the M model, belongs to the class of theories 
introduced by Seiberg [15] to give a dual description of conventional Af = 1 SQCD 
with the SU(iV c ) gauge group and Nf flavors of fundamental matter, where 

N c = N f -N 

(for reviews see Refs. [12, 16]). There are significant distinctions, however. 

Let us point out the main differences of the M model (3.5) from those introduced 
[15] by Seiberg: 

(i) Our theory has the U(iV) gauge group rather than SU(iV); 

(ii) Our theory has the FI D-term instead of a linear in M superpotential in 
Seiberg's models; 

(iii) We consider the case Nf = N which would correspond to Seiberg's iV c = 0. 
Our theory (3.5) is asymptotically free while Seiberg's dual theories give the most 
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reliable description of the original M = 1 SQCD in the range Nf < 3/2 N c which 
corresponds to Nf > 3N. In this range the theory (3.5) is not asymptotically free. 

In addition, it is worth noting that at Nf > N the vacuum (2.8), (2.10) be- 
comes metastable: supersymmetry is broken [17]. The N c = Nf — N supersymmetry- 
preserving vacua have vanishing VEV's of the quark fields and nonvanishing VEV 
of the M field 3 . The latter vacua are associated with the gluino condensation in 
pure SU(iV) theory, (AA) ^ 0, arising upon decoupling of Nf flavors [12]. In the case 
Nf = N considered here the vacuum (2.8), (2.10) preserves supersymmetry. Thus, 
despite a conceptual similarity between Seiberg's models and ours, dynamical details 
are radically different. 

To conclude this section let us mention that if a theory dual to the one in (3.5) 
were known our results would imply a non-Abelian confinement of quarks in the 
former theory. We will qualitatively discuss this issue in Sect. 8. 

4 Non-Abelian strings 

Non-Abelian strings were shown to emerge at weak coupling in Af = 2 supersymmetric 
gauge theories [5, 6, 7, 8]. The main feature of the non-Abelian strings is the presence 
of orientational zero modes associated with rotations of their color flux in the non- 
Abelian gauge group. This feature makes such strings genuinely non-Abelian. 

As long as the solution for the non-Abelian string suggested and discussed in 
[6, 7] for M = 2 SQCD does not depend on the adjoint fields it can be generalized to 
M = 1 SQCD upon introducing the mass term (2.2) for the adjoint fields and then 
taking the limit /ii j2 — > oo. This was done in Ref. [9]. However, as we have already 
explained above, M = 1 SQCD has the Higgs branch which obscures physics of the 
non-Abelian strings. The string becomes infinitely thick in the limit /Xj — > oo due to 
the presence of massless fields in the bulk. 

In particular, in [9] it turned out impossible to follow the fate of the confined 
monopoles (present in M = 2 SQCD) all the way down to M = 1 SQCD which one 
recovers in the limit = oo. Below we will show that this obstacle does not arise 
in the model (2.4). The reason is that M = 1 SQCD with the mesonic field M has 
no massless states in the bulk in the limit /ij — > oo, as was demonstrated in Sect. 3. 

Let us generalize the string solutions found in [6, 7] to the model (2.4). In 
addition to the conventional Abrikosov-Nielsen-Olesen (ANO) strings [3] this model 
supports Zn strings. These arise due to a nontrivial homotopy group 



3 This is correct for the version of the theory with ^-parameter introduced via superpotential. 




(4.1) 
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It is easy to see that this nontrivial topology amounts to winding of just one element 
of the diagonal matrix (2.8) at infinity. Such strings can be called elementary; their 
tension is 1/JV of that of the ANO string. The ANO string can be viewed as a bound 
state of N elementary strings. 

More concretely, the string solution (a progenitor of the non-Abelian string) 
can be written [6] as follows: 

/ <h(r) ... \ 

Q = , q = 0, 

... 2 (r) 

v ... e ia 0i(r) j 



A- 



SU(iV) 



/ 1 ... \ 



N 











V o o ... -(iv -l) J 



(dia) [-l + f NA (r)\, 



4 UW = ^ = Jf (dia)[l-f(r)], a = a a = M° = M a = 0, (4.2) 

where % — 1,2 labels coordinates in the plane orthogonal to the string axis, r and a 
are the polar coordinates in this plane and / is the unit N x N matrix. The profile 
functions <f>i(r) and 02 ( r ) determine the profiles of the scalar fields, while fNA{ r ) and 
f(r) determine the SU(A^) and U(l) fields of the string solution, respectively. These 
functions satisfy the following rather obvious boundary conditions: 

•Mo) = o, 

f NA (0) = l, /(0) = 1, (4.3) 

at r = 0, and 

M°°) = 02(00) = ^, 

f NA (oo) = 0, /(oo) = (4.4) 

at r = 00. 

As long as our ansatz (4.2) does not involve the fields q, a and M the classical 

string solution does not depend on Af = 2 SUSY breaking parameters. The classical 
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solution is the same as that found [6] in the M = 2 SQCD limit. In particular, the 
profile functions satisfy the following first-order equations: 

T Tr 0l(r) " Jf (/(r) + {N ~ 1)/jvA(r)) 0l(r) = ° ' 

r tr Mr) ~ IV (/(r) ~ fNA(r)) Mr) = ° ' 

~\ Tr m + [ (0l(r))2 + {N - 1} (02(r))2 - N *\ = ° ' 

^/^VA(r) + | [(<Mr)) 2 - (0 2 (r)) 2 ] = . (4.5) 
r dr 2 1 J 

Numerical solutions of the Bogomolny equations (4.5) for iV = 2 (Z 2 strings) were 
found in Ref. [6]. 

The string (4.8) is 1/2-BPS saturated. It automatically preserves two super- 
charges out of four present in the bulk theory. The tension of this elementary string 
is 

Ti = 27r£, (4.6) 
to be compared with the ANO string tension, 

T AN o = 2iV7r£ (4.7) 

in our normalization. 

The elementary strings are bona fide non-Abelian. This means that, besides triv- 
ial translational moduli, they acquire moduli corresponding to spontaneous break- 
ing of a non-Abelian symmetry. Indeed, while the "flat" vacuum (2.8), (2.9) and 
(2.10) is S\J(N)c+f symmetric, the solution (4.2) breaks this symmetry down to 
U(l)xSU(iV- 1). 

To obtain the non-Abelian string solution from the Zn string (4.2) we apply the 
diagonal color-flavor rotation (2.11) which preserves the vacuum. To this end it is 
convenient to pass to the singular gauge where the scalar fields have no winding at 
infinity, while the string flux comes from the vicinity of the origin. In this gauge we 
have (for details see the review paper [18]) 

q = 1[(JV -1)0 2 + + (0!-0 2 ) (n-n*- 1 
Af (N) = (n-n*-j^ £ij ^ f NA (r) , 
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A U(1) 



X 



0, a = a a : 



M° = M a = , 



(4.8) 



where we parametrize the matrices U of SU(iV) c+F rotations as follows: 

i 



1 

TV 



/ 1 



U 











\ 



u- 1 



-n l n* p + 



N p 



(4.9) 



... 1 
V ... -(N-l) J 

Here n l is a complex vector in the fundamental representation of SU(iV), and 

n\n l = \, (4.10) 

(l,p = 1,...,N are color indices). In Eq. (4.8) for brevity we suppress all SU(iV) 
indices. At n = {0, 1} we get the field configuration quoted in Eq. (4.2). 

The vector n l parametrizes orientational zero modes of the string associated with 
flux rotations in SU(iV). The presence of these modes makes the string genuinely non- 
Abelian. 

To derive an effective world-sheet theory for the orientational collective coordi- 
nates n l of the non-Abelian string we follow Refs. [6, 7, 19], see also the review [18]. 
From the string solution (4.2) it is quite clear that not each element of the matrix U 
will give rise to a modulus. The SU(iV— l)xU(l) subgroup remains unbroken by the 
string solution under consideration; therefore the moduli space is 



SU(iV) 



SU(JV- 1) x U(l) 



CP(N-l) 



(4.11) 



Assume that the orientational collective coordinates n l are slowly varying functions of 
the string world-sheet coordinates Xk, k — 0, 3. Then the moduli n l become fields of a 
(l+l)-dimensional sigma model on the world sheet. Since the vector n l parametrizes 
the string zero modes, there is no potential term in this sigma model. 

To obtain the kinetic term we substitute our solution, which depends on the 
moduli n l , in the action (2.4), assuming that the fields acquire a dependence on the 
coordinates x k via n\x k ). Then we arrive at the CP(N — 1) sigma model (for details 
see the review paper [18]), 



'CP(iV-l) 



2(3 J dtdz {(d k n*d k n) + (n*d k n) 2 } , (4.12) 
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where the coupling constant (3 is given by a normalizing integral defined in terms 
of the string profile functions. Using the first-order equations for the string profile 
functions (4.5) one can see that this integral reduces to a total derivative and given 
by the flux of the string determined by /jva(O) = 1, namely 

= 5- (4-13) 

The two-dimensional coupling constant is determined by the four- dimensional non- 
Abelian coupling. 

The relation between the four- dimensional and two-dimensional coupling con- 
stants (4.13) is obtained at the classical level. In quantum theory both couplings 
run. So we have to specify a scale at which the relation (4.13) takes place. The 
two-dimensional CP(N — 1) model is an effective low-energy theory good for the 
description of internal string dynamics at low energies, much lower than the inverse 
thickness of the string which, in turn, is given by (fev^- Thus, g^VC plays the role of a 
physical ultraviolet cutoff in (4.12). This is the scale at which Eq. (4.13) holds. Below 
this scale, the coupling f3 runs according to its two-dimensional renormalization-group 
flow. 

The sigma model (4.12) is asymptotically free [20]; at large distances (low en- 
ergies) it gets into the strong coupling regime. The running coupling constant as a 
function of the energy scale E at one loop is given by 

4n[3 = Nln ( —— ]+•••, (4.14) 

\Acp(jv-i) / 

where Acp(n-i) is the dynamical scale of the CP(N — 1) model. As was mentioned 
above, the UV cut-off of the sigma model at hand is determined by g<i \J\. Hence, 

AS P{ N-i ) = 9?t N/2 e~ 1 *- (4-15) 

Note that in the bulk theory, due to the VEV's of the squark fields, the coupling 
constant is frozen at g2\f£,- There are no logarithms in the bulk theory below this 
scale. Below <?2"\/£ the logarithms of the world-sheet theory take over. 
In the limit of large Hi we are interested in here, 



A*2 > 92\[l, 



the coupling constant #2 of the bulk theory is determined by the scale Ajv"=i of the 
M = 1 SQCD (3.5) with the M field included, as shown in Eq. (3.6). In this limit 
Eq. (4.15) gives 

AcP(iv-i) = ^% , (4.16) 
92V C 
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where we take into account that the first coefficient of the (3 function in M = 1 SQCD 
is 2N. 

To conclude this section let us note a somewhat related development: non-BPS 
non-Abelian strings were recently considered in metastable vacua of a dual description 
of N = 1 SQCD at N f > N in Ref. [21]. 

5 Fermionic sector of the world-sheet theory 

In this section we discuss the fermionic sector of the low-energy effective theory on 
the world sheet of the non-Abelian string in Af = 1 SQCD with the M field, as well 
as supersymmetry of the world-sheet theory. First we note that our string is 1/2 BPS 
suturated. Therefore in the M = 2 limit ( when M = 2 breaking parameters /ij and 
h vanish) four supercharges out of eight present in the bulk theory are automatically 
preserved on the string world sheet. They become supercharges in the CP(N — 1) 
model (4.12). 

For simplicity in this section we will discuss the case N = 2 limiting ourselves 
to the CP(1) model. Generalization to arbitrary N is straightforward. The action of 
the (2,2) supersymmetric CP(1) model model is 

Scili) = pfdtdzfyd k S a ) 2 + ±Xii(do-id 3 )xi 

+ \x a 2^{^,+^^,)x a 2-\{x a l X a 2) 2 }, (5-1) 

where we used the fact that CP (I) is equivalent to the 0(3) sigma model defined in 
terms of a unit real vector S a , 

S a = n\T a n\ (S a ) 2 = 1 . (5.2) 

This model has two real bosonic degrees of freedom. Two real fermion fields Xi an d 
%2 are subject to constrains 

X\S a = 0, X a 2 S a = . (5.3) 

Altogether we have four real fermion fields in the model (5.1). 

Now we break M = 2 supersymmetry of the bulk model by switching on pa- 
rameters ^ and h. The l/2-"BPS-ness" of the string solution requires only two 
supercharges. However, as we will show below, the number of the fermion zero modes 
on the string does not change. This number is fixed by the index theorem. Thus, the 
number of (classically) massless fermion fields in the world sheet CP(N — 1) model 
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does not change. It was shown in [9] that the (2,2) supersymmetric sigma model 
with the CP(N — 1) target space does not admit (0,2) supersymmetric deformations. 
Therefore, it was concluded in [9] that the world sheet theory has "accidental" SUSY 
enhancement. A similar phenomenon was found earlier in [22] for domain walls. 

On the other hand, in the recent publication [23] it was suggested that super- 
orientational zero modes can mix with supertranslational modes. It was shown that 
the sigma model with the C x CP(N — 1) target space does admit (0,2) supersym- 
metric deformations. It is not clear at the moment if this mixing really occurs in the 
effective theory on the string. If it occurs then the emerging (0,2) supersymmetric 
C x CP(N — 1) model has a /i-deformed four-fermion interaction 



°CP{1) 



= pjdtdz^(d k S a ) 2 + ±Xii(d -id 3 )x a i 



+ lx a 2i(do + id 3 ) X a 2-l 



2 i + cMV(<fe 2 



(5.4) 



where c is an unknown coefficient. Also the first constraint in (5.3) is replaced with 
X\S a — c/2 (/X2C1 + A^Ci), where £i is the right-moving two-dimensional fermion field 
associated with the supertranslational zero modes. If this conjecture [23] is correct 
the four-fermion term disappears in the large-// limit. To find out which scenario is 
correct one has to calculate the coefficient in front of the four-fermion term in (5.4). 
We left this for future work. 

In any case, the world sheet supersymmetric model has N vacua which are 
interpreted as N elementary strings of the bulk theory. This number is protected by 
Witten's index and survives M = 2 breaking deformations. We will use this result 
in the next section. The kinks which interpolate between these vacua are confined 
monopoles. Below we will show that the occurrence of four (4(iV — 1) in the general 
case) superorientational fermion zero modes on the non-Abelian strings follows from 
an index theorem. In Appendix we present explicit solutions for these modes for the 
case N = 2. 



5.1 Index theorem 

The fermionic part of the action of the model (3.5) is 



S 



ferm 



d A x { -^\ a ip\ a + -3 A^A + Tr i>i 



92 



9i 



Tr 



2i r__ i % r ~ - ~ 

+ -Tr [0C\ + ^= Tr [g(A^) - (^A)g + (^A)g - q(X^) 
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+ ^Tr 



+ iTr 



q2T a (\ a *P) - {^\ a )2T a q + {^\ a ) 2T a q - q2T a (A» 

9(^0 + (HQ + mo + q($£) 



+ i Tr (ipipM + ipipM^j } , 



(5.5) 



where the matrix color-flavor notation is used for matter fermions {jp a ) kA and {jp a )Ak 
and the traces are performed over the color-flavor indices. Contraction of the spinor 
indices is assumed inside all parentheses, for instance, {Xiji) = \ a ip a ■ Moreover, ( 
denotes the fermion component of matrix M superfield, 



(5.6) 



In order to find the number of the fermion zero modes in the background of the 
non-Abelian string solution (4.8) we have to carry out the following program. Since 
our string solution depends only on two coordinates X{ (i — 1, 2), we can reduce our 
theory to two dimensions. Given the theory defined on the (xi,X2) plane we have to 
identify an axial current and derive the anomalous divergence for this current. In two 
dimensions the axial current anomaly takes the form 

dija ~ F*, (5.7) 

where F* = (l/2)e i jF i j is the dual U(l) field strength in two dimensions. 

Then integral over the left-hand side over the (xi, X2) plane gives us the index of 
the 2D Dirac operator v coinciding with the number of the 2D left-handed minus 2D 
right-handed zero modes of this operator in the given background field. The integral 
over the right-hand side is proportional to the string flux. This will fix the number 
of the chiral fermion zero modes 4 of the string with the given flux. Note that the 
reduction of the theory to two dimensions is an important step in this program. The 
anomaly relation in four dimensions involves the instanton charge F*F rather than 
the string flux and is therefore useless for our purposes. 

The reduction of Af = 1 gauge theories to two dimensions is discussed in detail 
in [24] and here we will be brief. Following [24] we use the rules 

r ^or,^ + ), r->(^-,^ + ), 

A*^(A-,A + ), C a -(C,C + )- (5.8) 

4 Chirality is understood as the two-dimensional chirality. 
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Field 


^+ 


if). 


^+ 


^_ 


A+ 


A_ 


C+ 


C- 


q 


q 


U(l)# charge 


-1 


1 


-1 


1 


-1 


1 


-1 


1 








U(l)^j charge 


-1 


1 


1 


-1 


-1 


1 


1 


-1 









Table 1: The U(l)i? and U(l)^ charges of fields of the two-dimensional reduction of the theory. 

With these rules the Yukawa interactions in (5.5) take the form 

^Yukawa = iV2Tr -q{\_^ + - \ + if>_) + (t/>_A + - ^+A_)g + c.c. 

- i Tr [q(ip-( + - V'+C-) + $-qC+ ~ + cc] , (5.9) 

where the color matrix A = (1/2) A + T a \ a . 

It is easy to see that ^Yukawa is classically invariant under the chiral U(1)r 
transformations with the U(l)^ charges presented in Table 1. The axial current 
associated with this U(l)# is not anomalous [24]. This is easy to understand. In 
two dimensions the chiral anomaly comes from the diagram shown in Fig. 1. The 
U(l)# chiral charges of the fields if) and if) are the same while their electric charges 
are opposite. This leads to cancellation of their contributions to this diagram. 

It turns out that for the particular string solution we are interested in the classi- 
cal two-dimensional action has more symmetries than generically, for a general back- 
ground. To see this, please, note that the field q vanishes on the string solution (4.8). 
Then the Yukawa interactions (5.9) reduce to 



iV2Tr \-q(\^+ ~ K4>- 



i Tr 



- i>+q(- 



+ c.c. 



(5.10) 



The fermion if) interacts only with A's while the fermion if) interacts only with (. Note 
also that the interaction in the last line in (5.5) is absent because M = on the string 
solution. This property allows us to introduce another chiral symmetry in the theory, 
the one which is relevant for the string solution. We will refer to this extra chiral 
symmetry as U(l)^. 

The U(l)^, charges of our set of fields are also shown in Table 1. Note that if> 
and if) have the opposite charges under this symmetry. The corresponding current 
then has the form 



Jib 



\ 



if)_if)_ - if)+if)+ - if)_if)- + H 

V —iif)-if)- — iif)+if)+ + iif)_if)- + iif) + if) + + ■ ■ ■ j 



(5.11) 
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Figure 1: Diagram for the chiral anomaly in two dimensions. The solid lines denote fermions tp, 
if), the dashed line denotes photon, while the cross denotes insertion of the axial current. 

where the ellipses stand for terms associated with the A and ( fields which do not 
contribute to the anomaly relation. 

Clearly, in quantum theory this symmetry is anomalous. Now the contributions 
of the fermions ip and ip double in the diagram in Fig. 1 rather than cancel. It is not 
difficult to find the coefficient in the anomaly formula 

N 2 

dd i5 = —F*, (5.12) 

71 

which can be normalized e.g. from [25]. The factor N 2 appears due to the presence 
of 2N 2 fermions ip k and ipAk- 

Now, taking into account that the flux of the Zn string under consideration is 

Jd 2 xF* = ^, (5.13) 

(see the expression for the U(l) gauge field for the solution (4.2) or (4.8)) we conclude 
that the total number of the fermion zero modes in the string background 

v = AN . (5.14) 

This number can be decomposed as 

v = AN = A(N- 1) +4, (5.15) 

where 4 is the number of the supertranslational modes while 4(7V — 1) is the number 
of the superorientational modes. Four supertranslational modes are associated with 
four fermion fields in the two-dimensional effective theory on the string world sheet, 
which are superpartners of the bosonic translational moduli Xq and Uq. Furthermore, 
4(iV — 1) corresponds to 4(iV — 1) fermion fields in the M = 2 CP(N — 1) model 
on the string world sheet (4.12). CP(N — 1) describes dynamics of the orientational 
moduli of the string. For N = 2 the latter number (A(N — 1) = 4) counts four fermion 
fields Xii X2 m the model (5.1) or (5.4). 

We explicitly determine four superorientational fermion zero modes for the case 
N = 2 in Appendix. Note that the fermion zero modes of the string in M = 1 SQCD 
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with the M field are perfectly normalizable provided we keep the coupling constant 
h nonvanishing. Instead, in conventional M = 1 SQCD without the M field the 
second pair of the fermion zero modes (proportional to xl ) become non-normalizable 
[9]. This is related to the presence of the Higgs branch and massless bulk states in 
conventional M = 1 SQCD. As was already mentioned more than once, in the M 
model, Eq. (3.5), we have no massless states in the bulk. 

Note that in both translational and orientational sectors the number of the 
fermion zero modes is twice larger than the one dictated by l/2-"BPS-ness." 

6 Evolution of the monopoles 

Since supersymmetric CP(N — 1) model is an effective low-energy theory describing 
world sheet physics of the non-Abelian string, all consequences of this model ensue, 
in particular, N degenerate vacua and kinks which interpolate between them — the 
same kinks that we had discovered in M = 2 SQCD [7] and interpreted as (confined) 
non-Abelian monopoles, the descendants of the 't Hooft-Polyakov monopole [26]. 

Let us briefly review the reason for this interpretation [27, 7, 8] and discuss what 
happens with these monopoles as we deform our theory and eventually end up with 
the M model. It is convenient to split this deformation into several distinct stages. 
We will describe what happens to the monopoles as one passes from one stage to 
another. 

A qualitative evolution of the monopoles under consideration as a function of 
the relevant parameters is presented in Fig. 2. 

• We start from Af = 2 SQCD turning off the M = 2 breaking parameters h 
and /i's as well as the FI parameter in the theory (2.4), i.e. we start from the 
Coulomb branch of the theory, 

/i! = /i 2 = 0, h = 0, f = 0, M^O. (6.1) 

As was explained in Sect. 2, the field M is frozen in this limit and can take 
arbitrary values (the notorious flat direction) . The matrix plays the role of 
fixed mass parameter matrix for the quark fields. First we consider the diagonal 
matrix M, with distinct diagonal entries, 

M^ = diag{M 1 ,...,M^}. (6.2) 

Shifting the field a one can always make J2a — in the limit [i\ = 0, 
therefore M° = 0. If all M^'s are different the gauge group SU(iV) is broken 
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1=0 
Am^O 




Confined monopole, 
quasiclassical regime 



Confined monopole, 
highly quantum regime 



Figure 2: Various regimes for the monopoles and flux tubes in the simplest case of two flavors. 



down to U(1)( 7V " 1 ) by a VEV of the SU(7V) adjoint scalar 

(a?) = --^=5?M l . 



(6.3) 



Thus, there are 't Hooft-Polyakov monopoles embedded in the broken gauge 
SU(iV). Classically, on the Coulomb branch the masses of (N — 1) elementary 
monopoles are proportional to 

\{M A -M A+1 )\/g 2 2 

This is shown in the upper left corner of Fig. 2 for the case 

N = 2, Am = Mx-M 2 . 

In the limit (Ma — Ma+x) — * the monopoles tend to become massless, for- 
mally, in the classical approximation. Simultaneously their size become infinite 
[28]. The mass and size are stabilized by confinement effects which are highly 
quantum. The confinement of monopoles occurs in the Higgs phase, at £ ^ 0. 

Now we introduce the FI parameter £ which triggers the squark condensation. 
The theory is in the Higgs phase. We still keep M = 2 breaking parameters h 
and /i's vanishing, 



fil = fi 2 = Q 1 h = 0, £ ^ 0, M ^ 0. 



(6.4) 
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If we allow £ to be nonvanishing but small, 

|M A |»^, (6.5) 

then the effect which comes into play first is the spontaneous breaking of the 
gauge SU(iV) by the condensation of the adjoint scalars. The next gauge sym- 
metry breaking, due to £ 7^ 0, which leads to complete Higgsing of the model 
and the string formation (confinement of monopoles) is much weaker. Thus, we 
deal here with the formation of "almost" 't Hooft-Polyakov monopoles, with a 
typical size ~ \{M A — M A+1 ) | _1 . Only at much larger distances, ~ £ -1//2 , the 
chromoelectric charge condensation enters the game, and forces the magnetic 
flux, rather than spreading evenly a la Coulomb, to form flux tubes (the upper 
right corner of Fig. 2). 

Let us verify that the confined monopole is a junction of two strings. At M A 7^ 
the global SV(N) c +f group is broken by condensation of the adjoint scalars 
(6.3), and non-Abelian strings become Abelian strings. Their orientational 
moduli space is lifted [7, 8]. Consider the junction of two Zn strings (4.8), 
namely A-tYi string with 

n l = 5 l A (6.6) 
and "neighboring" (A + l)-th string with 

n l = 5 l A+1 , (6.7) 

(cf. solution (4.2) which is written for A + 1 = N). The flux of the junction 
is given by the difference of the fluxes of these two strings. Using (4.8) we get 
that the flux of the junction is 

4tt x diagi {...0, 1, -1, 0, ...} (6.8) 

with the nonvanishing entries located at the positions A and (A + 1). These 
are exactly the fluxes of N — 1 distinct 't Hooft-Polyakov monopoles occurring 
in the SU(iV) gauge theory provided that SU(iV) is spontaneously broken down 
to L^l)^ -1 . We see that in the quasiclassical limit of large \M A \ the Abelian 
monopoles and the junctions of the Abelian Zn strings are in one-to-one corre- 
spondence. 

At large M A the monopoles, albeit confined, are weakly confined. Now, if 
we further reduce \M A \, 

Acp(n-i) < \M A \ < y^, (6.9) 
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the size of the monopole ~ \{M A — M A+1 ) | _1 becomes larger than the trans- 
verse size of the attached strings. The monopole gets squeezed in earnest by 
the strings — it becomes a bona fide confined monopole (the lower left corner 
of Fig. 2). At nonzero M A the effective CP(N — 1) model on the string world 
sheet becomes massive with the potential determined by so called twisted mass 
terms [27, 7, 8]. Two Z N strings corresponds to two "neighboring" vacua of the 
CP(N — 1) model . The monopole (aka the string junction of two Z N strings) 
is interpreted as a kink interpolating between these two vacua. 

In [7] the first order equations for the 1/4 BPS string junction of two Z 2 
strings were explicitly solved in the case N = 2, and the solution shown to 
correspond to the kink solution of the two-dimensional CP(1) model. Moreover, 
it was shown that the mass of the monopole matches the mass of the CP(1)- 
model kink both in the quasiclassical (Am A C p(i)) and quantum (Am <C 
Acp(i)) limits. 

• Now let us switch off the mass differences M A still keeping the M = 2 breaking 
parameters vanishing, 

^=^2 = 0, h = 0, f^O, M = 0. (6.10) 

The values of the twisted mass in CP(N — 1) model equal M A while the size 
of the twisted-mass sigma- model kink/confined monopole is of the order of 
~\(M A -M A+1 )\- 1 . 

As we further diminish M A approaching A C p(N-i) an d then getting below 
A-cp(n-i), the size of the monopole grows, and, classically, it would explode. 
This is where quantum effects in the world-sheet theory take over. It is natural 
to refer to this domain of parameters as the "regime of highly quantum dynam- 
ics." While the thickness of the string (in the transverse direction) is ~ C" 1 ^ 2 , 
the z-direction size of the kink representing the confined monopole in the highly 
quantum regime is much larger, ~ A-cp(jv-i)> see the lower right corner in Fig. 2. 

In this regime the confined monopoles become truly non-Abelian. They no 
longer carry average magnetic flux since 

(n') = 0, (6.11) 

in the strong coupling limit of the CP(N — 1) model [30]. The kink/monopole 
belongs to the fundamental representation of the S\J(N) c +f group [30, 31]. 

Let us stress that in the limit M A = the global group S\J(N) C+F is restored 
in the bulk and both strings and confined monopoles become non-Abelian. One 



23 



might argue that this restoration could happen only at the classical level. One 
could suspect that in quantum theory a "dynamical Abelization" ( i.e. a cascade 
breaking of the gauge symmetry U(iV)— >U(l) iV — > descrete subgroup ) might 
occur. This could have happened if the adjoint VEV's that are classically vanish 
at M = (see (2.9)) could have developed dynamically in quantum theory. 

At Ma 7^ the global SXJ(N)c+f group is explicitly broken down to L^l)^ -1 
by quark masses. At Ma = this group is classically restored. If still it could 
be dynamically broken this would mean a spontaneous symmetry breaking. 

Let us show that this does not happen in the theory at hand. First of all, if a 
global symmetry is not spontaneously broken at the tree level then it cannot 
be broken by quantum effects at week coupling in "isolated" vacua. Second, if 
the global group S\J(N) c +f were broken spontaneously at Ma = this would 
ensure the presence of massless Goldstone bosons. However, we know that there 
are no massless states in the spectrum of the bulk theory, see Sect. 2, 3. 

Finally, the breaking of SU(A^) C+ ^ in the Ma = limit would mean that 
the twisted masses of the world sheet CP(N — 1) model would not be given 
by Ma', instead they would be shifted, = M a + ca^-cp(n-i), where ca 

are some coefficients. In [7, 8] it was shown that the BPS spectrum of the 
CP(N — 1) model on the string should coincide with the BPS spectrum of 
the four- dimensional bulk theory on the Coulomb branch because the central 
charges which determine masses of the BPS states cannot depend on the non- 
holomorphic parameter £. The BPS spectrum of the CP(N — 1) model is 
determined by while the BPS spectrum of the bulk theory on the Coulomb 
branch is determined by Ma- In [32] it was shown that the BPS spectrum of 
both theories coincide at = Ma- Thus, we conclude that ca = and the 

twisted masses go to zero in the Ma = limit. Again we conclude that the 
global SXJ(N)c+f group is not broken in the bulk and both strings and confined 
monopoles become non-Abelian at Ma = 0. 

• Thus, at zero Ma we still have confined "monopoles" (interpreted as kinks) 
stabilized by quantum effects in the world-sheet CP(N — 1) model. Now we 
can finally switch on the Af = 2 breaking parameters //j and h, 

Hi^O, h^O, ^0, M = 0. (6.12) 

Note that the last equality here is automatically satisfied in the vacuum, see 
Eq. (2.10). 

As we discussed in Sects. 4 and 5 the effective world-sheet description of 
the non-Abelian string is still given by supersymmetric CP(N — 1) model. This 
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model obviously still has iV vacua which should be interpreted as N elementary 
non-Abelian strings in the quantum regime, and BPS kinks can interpolate 
between these vacua. These kinks should still be interpreted as non-Abelian 
confined monopoles/string junctions. 

Note that although the adjoint fields are still present in the theory (2.4) their 
VEV's vanish (see (2.9)) and the monopoles cannot be seen in the semiclassical 
approximation. They are seen as the CP(N — 1) model kinks. Their mass and 
inverse size is determined by A.qp(n-i) which in the limit of large //j is given by 
Eq. (4.16). 

• Now, at the last stage, we take the limit of large masses of the adjoint fields in 
order to eliminate them from the physical spectrum altogether, 

Atj^oo, h^O, f^O, M = 0. (6.13) 

The theory flows to M = 1 SQCD extended by the M field. 

In this limit we get a remarkable result: although the adjoint fields are elim- 
inated from our theory and the monopoles cannot be seen in any semiclassical 
description, our analysis shows that confined non-Abelian monopoles still exist 
in the theory (3.5). They are seen as CP(N — l)-model kinks in the effective 
world-sheet theory on the non-Abelian string. 

7 A brane perspective 

Let us elucidate how some important features of the consideration above are seen in 
the brane picture. To this end we will rely on Type IIA string approach to our M 
model. Consider the brane picture for Af = 2 and Af = 1 SQCD (see Ref. [33] for 
a review). We will limit ourselves to the special case of equal numbers of colors and 
flavors relevant to the present work. 

The Af = 2 theory involves N D4 branes extended in the directions of the (0,1, 
2, 3, 6) coordinates, two NS5 branes with coordinates along (0, 1, 2, 3, 4, 5), localized 
at Xq = and xq = 1/g 2 and Nf = N D6 branes with the world volume along (0, 1, 
2, 3, 7, 8, 9). The D4 branes are stretched between NS5 branes along xq, while the 
coordinates of D6 branes in xq are arbitrary. The NS5 branes can be split in the x^ 
direction which corresponds to the introduction of the Fayet-Iliopoulos term in the 
U(l) factor of U(iV), namely, 

8x 7 = £. 

The Higgs branch in this theory occurs when the D6 branes touch the D4 branes. 
After this, the D4 branes can split in pieces which can be moved in the (7, 8, 9) 
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directions. The coordinates of these pieces in in the (7, 8, 9) directions, along with 
the Wilson line of A 6 , yield coordinates on the Higgs branch of the moduli space. 
Fluctuations of the ends of the D4 branes in the (4, 5) plane provide the coordinates 
on the Coulomb branch of the moduli space. 

To break M = 2 SUSY down to M = 1 we rotate one of the NS5 branes. The 
angle of rotation corresponds to the mass of the adjoint scalar in the superpotential 
(2.2). The fact that this superpotential does not vanish removes the Coulomb branch 
of the moduli space. The positions of the D4 branes in the (4,5) plane are now fixed. 

Now, let us switch on the meson field M. It turns out that it emerges as a 
particular limiting brane configuration in the setup described above, without any 
additional branes. Consider the situation when the Xq coordinates of all D6 branes 
are the same. First, in this limit the open strings connecting the pairs of the D6 branes 
yield a massless field which is in the adjoint representation with respect to the flavor 
group U(iV). In the field-theory language this is nothing but our M field. Taking into 
account the standard three-string vertex we immediately derive the superpotential 
Wm = Tr M Q Q. On the other hand, since all D6 branes have the same x 6 coordinate, 
it is impossible to split the pieces of the D4 branes — such a splitting would require 
different values of xq for the pair of the D6 branes. Thus, the Higgs branch disappears. 
We see that in the brane language the introduction of the M field is in one-to-one 
correspondence with the disappearance of the Higgs branch. 

Consider now the evolution of the monopoles discussed in Sect. 6 within the 
framework of the brane picture. In the M = 2 theory in the regime (6.1) the monopole 
is represented by a D2 brane stretched between two NS5 branes in the x 6 direction 
and two D4 branes located at x^a = Ma and x^a+i) = M^+i), which yields the 
correct monopole mass 

\(M A -M A+1 )\ 

gl 

Switching on the Fayet-Iliopoulos parameter parameter £ in the regime (6.4) 
corresponds to a displacement of one of the NS5 branes in the xj direction. Since the 
D4 branes split in two pieces at the common x G coordinate where the D6 branes are 
located, and each piece is attached to the NS5 brane, a squark condensate develops. 
It is proportional to This regime supports quasiclassical non-Abelian strings 
which have a transparent geometrical interpretation [5] . The non-Abelian strings are 
identified with the D2 brane parallel to the D6 branes stretched between two NS5 
branes along the x 7 coordinate. Geometrically, the string tension equals 5x 7 , in full 
agreement with the field-theory result. 

The D2 brane representing the monopole in the Higgs phase is located as follows. 
It extends along two coordinates, x 6 and x 4 . Along the x 6 coordinate the D2 brane 
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is stretched between the common position of the D6 branes and the NS5 brane. In 
the £4 direction it is stretched between two D4 branes. 

From this picture one immediately recognizes the monopole to be a junction of 
two non-Abelian strings since it is stretched between two different non-Abelian strings 
in the £4 direction. If one switches off the Fayet-Iliopoulos term then the monopole 
in the Higgs phase geometrically smoothly transforms into the 't Hooft-Polyakov 
monopole. 

This picture implies that in the semiclassical regime of large Ma the monopole 
mass is the same as the mass of the 't Hooft-Polyakov monopole. With Ma decreasing 
we eventually find ourselves in the purely quantum regime described by lifting type 
IIA string to M-theory and, hence, lifting the D2 brane to M2 brane. In M-theory 
the monopole in the Higgs phase can be easily described by the M2 brane wrapping 
the appropriate circle on the Riemann surface, using its identification with the kink 
in CP(JV- 1) model [32]. 

Finally in the regime (6.12) we rotate one of the NS5 branes which results in 
vanishing vacuum expectation values of the adjoint scalars. However, the M2 brane 
representing the non-Abelian string is still clearly identified. The monopoles are the 
M2 branes wrapped around the Riemann surface responsible for this regime upon 
rotation of the branes. 

Let us emphasize that the monopoles in all regimes are represented by the M2 
branes, and their evolution from the Coulomb branch to the Higgs one corresponds 
just to different placement of one and the same brane in a certain brane background. 

Note that the brane picture suggests the possibility of a more general situation, 
when only k of the D6 branes have the same xq coordinates. Then, the massless 
meson field M belongs to the U(/c) subgroup of the flavor group. In particular, one 
can consider the case Nf > N, introduce a meson field of some rank and perform the 
standard Seiberg duality transformation by exchanging two NS5 branes. 

8 Discussion and conclusions 

Let us summarize our findings. Deformation of M = 2 SQCD leads us to the M 
model, M = 1 SQCD supplemented by the M field, see (3.5). We observe confined 
non-Abelian monopoles in this model which has no monopoles whatsoever in the 
semiclassical limit. Why we are sure that the objects we observe are "non-Abelian 
monopoles"? We know this because we can start from the Af = 2 theory on the 
Coulomb branch were the standard 't Hooft-Polyakov monopoles are abundant, and 
trace their evolution stage by stage, as one varies the adjustable parameters to even- 
tually arrive at M = 1 SQCD. 
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This is the main result of the present paper. As was mentioned above the 
confined monopoles are in the highly quantum regime so they do not carry average 
magnetic flux (see Eq. (6.11)). They are genuinely non-Abelian. Moreover, they 
acquire global flavor quantum numbers. In fact, they belong to the fundamental rep- 
resentation of the global SU(iV)c + .F group (see Refs. [30, 31] where this phenomenon 
is discussed in the context of the CP(N — l)-model kinks). 

In particular, the monopole-antimonopole "meson" formed by the string config- 
uration shown in Fig. 3 belongs to the adjoint representation of the global "flavor" 
group SU(iV)c + .F, in accordance with our expectations. Similar there are "baryons" 
built of N monopoles connected by strings to each other to form a close necklace 
configuration. 




Figure 3: Monopole and antimonopole bound by strings in a meson. Open and closed circles 
denote monopole and antimonopole, respectively. 

We believe that the emergence of these non-Abelian monopoles can shed light 
on mysterious objects introduced by Seiberg: "dual magnetic" quarks which play an 
important role in the description of N = 1 SQCD at strong coupling [15, 12]. 

It is curious to note that monopole-like configurations apparrently occur in lattice 
QCD. In particular, in the recent publications [34] the occurence of the monopole- 
like configurations is traced back to the color-octet operator q~T a q. We would like to 
stress that the non-Abelian monopoles observed here are totally different. In the limit 
fi — > oo all traces of "Abelization" ( i.e. cascade breaking of the gauge symmetry 
U(N)— >U(1) 7V — > descrete subgroup ) typical of the M = 2 limit are erased! In fact, 
it is clear from (2.8) that (q~T a q) = in the M-model vacuum and cannot be used 
to construct monopoles. Our monopoles are not seen classically. The confined non- 
Abelian monopoles emerge as CP(N — l)-model kinks living on the string, deep in 
the quantum regime. 

Now, let our imagination run away with the hypothetical dual of the M model. 
In this model it is not chromomagnetic but rather chromoelectric flux tubes that 
will form (upon "monopole" condensation) in a highly quantum regime. The number 
of degenerate chromoelectric flux tubes must grow with N. Quarks are confined; 
inside mesons a quark and its anti-partner must be attached to a pair of strings, in 
contradistinction with QCD where the confining bond between quark and anti-quark 
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is built from a single string. It is thus clear that even if a dual to the M model is 
found, it is not yet quite QCD. However, it is pretty close. 
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Appendix: Superorientational zero modes 

In this Appendix we find explicit expressions for four superorientational fermion zero 
modes of the non-Abelian string in the theory (3.5) with N = 2. Half-criticality of the 
string in question ensures that two supercharges are preserved in the world-sheet the- 
ory. Following the general method of [7, 9] we generate two superorientational fermion 
zero modes applying SUSY transformations to our string solution (4.8). Essentially 
repeating the calculation made in [9] we get 



V I J Ak 20 2 

sal _ * X\-iX2 01 \ a . a bc qb ,.c' 

A - /jt o J NA T X2+ l£ b X2 > 

V2 r 2 <p2 1 J 

\ a2 = 0. (A.l) 

We see that supersymmetry generates for us only two fermion superorientational 
modes out of four predicted by the index theorem. They are parametrized by the two- 
dimensional fermion field \%- This was expected, of course. The modes proportional 
to Xi do not appear. 

The nonzero fermion fields here have the U(l)^ charge +1 while the fields which 
are zero have charge —1. Clearly we need to find two more zero modes of charge 
+ 1. We do it by explicitly solving the Dirac equations. From the fermion part of the 
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action of the model (5.5) we get the relevant Dirac equations 

^A a + 4fTr^r a g = 0, 

l^ + l T r^r a = 0, 

t^+ i -C(qr a ) = 0. (A.2) 

After some algebra one can check that equations (A.l) do satisfy the first and the 
third of the Dirac equations (A.2) provided the first-order equations for the string 
profile functions (4.5) are satisfied. 

Now let us find two additional fermion zero modes by solving the second and 
the fourth of the Dirac equations (A.2). The fields with the U(l)^ chiral charge — 1 
vanish, 

^ =0, C al = 0. (A.3) 

In order to find the fields with the U(l) ^ chiral charge +1 we use the following ansatz 5 
(cf. Ref. [9]), 



= X ~^m [Xt + ie^xt]. (A.4) 

Here we introduce two profile functions ((r) and ip(r) parameterizing the fermion 
fields ( 2 and tpi- 

Substituting (A.4) into the Dirac equations (A.2) we get the following equations 
for fermion profile functions: 

^ + -i> - ^(/ + Ina)^ + % 0i C = o, 

ar r 2r 

- i^+ i ^^ = °- ( A - 5 ) 

Below we present the solution to these equations in the limit 

h^gf-gl (A.6) 



5 One can show that profile functions in front of all other possible structures have singular behavior 
either at r — or at r = oo. 
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This latter assumption is not a matter of principle, rather it is just a technical point. 
It allows us to find an approximate analytic solution to Eqs. (A. 5). If the condition 
(A. 6) is met the mass of the fermions if) and (, 

(A.7) 

(see Eqs. (3.2) and (3.3)) becomes much smaller than the masses of the gauge bosons 
(see Eqs. (2.13) and (2.14); note that the fermionsj0 and ( are superpartners of q 
and M and have the same mass). Thus, the fields if) and ( develop long range tails 
with the exponential fall-off determined by the small masses (A.7). This allows us to 
solve Eqs. (A. 5) analytically treating separately the domains of large and small r. 

In the large r domain, at r ^> mw, we can drop the terms in (A. 5) containing / 
and /na and use the first equation to express if) in terms of (. We then get 

2i (I t . A . 

Substituting this into the second equation in (A. 5) we obtain 

^C + ^C-^C = 0. (A.9) 
dr z r ar 

This is a well-known equation for a free field with mass mo in the radial coordinates. 
Its solution is well-known too, 

( = - i ±^K (m r), (A.10) 

where K (x) is the imaginary-argument Bessel function, and we fix a certain con- 
venient normalization (in fact, the normalization constant of the profile functions is 
included in xD- At infinity K Q (x) falls-off exponentially, 

K (x)~^=, (A.ll) 

while at x — > it has a logarithmic behavior, 

K (x)~hi-. (A. 12) 

x 

Taking into account Eq. (A. 8) we get the solutions for the fermion profile functions 

at r > l/m w , 

( = - l ±^K (m r), ^ = -± Ko (m r). (A.13) 
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In particular, at r <C l/m we have 

C~-^ln— , (A.14) 

In the intermediate domain r < l/m w we neglect the small mass terms in (A. 5). 
We then arrive at 

^ + ^-^(/ + /^W = 0. (A.15) 

The first equation here shows that ( =const, while the second one is identical to the 
equation for the string profile function 1; see Eq. (4.5). This gives the fermion profile 
functions at intermediate r, 

C = 4^m^, = -L, (A.i6) 

where we fixed the normalization constants matching this solutions with the ones in 
the large-r region, see (A.14). 

Equations (A. 13) and (A. 16) present our final result for the fermion profile func- 
tions. They determine two extra fermion superorientational zero modes proportional 
to Xi via Eq. (A.4). 

Now if we substitute the fermion zero modes (A.l) and (A.4) in the action (5.5) 
we get the effective Af = 2 CP(1) model (5.1) on the string world sheet, 6 cf. Ref. [9]. 
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